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Abstract. This paper is concerned with the bifurcation of limit cycles in 
general quadratic perturbations of quadratic codimcnsion-four centers Q4. 
Gavrilov and Iliev set an upper bound of eight for the number of limit cycles 
produced from the period annulus around the center. Based on Gavrilov-Iliev's 
proof, we prove in this paper that the perturbed system has at most five limit 
cycles which emerge from the period annulus around the center. We also show 
that there exists a perturbed system with three limit cycles produced by the 
period annulus of Q4. 



1. Introduction and statement of the main results 

In this paper we study the bifurcation of limit cycles in plane quadratic systems 
under small quadratic perturbations. We assume that the unperturbed system 
has at least one center. Taking a complex coordinate z — x + iy and using the 
terminology from [25], the list of quadratic centers at z = looks as follows: 



2 -J- 2\z\ 2 + (b + ic)z 2 , Hamiltonian (Qf), 



z = —iz — z 

z = —iz + az 2 + 2\z\ 2 + bz 2 , reversible {Qf), 

z = —iz + Az 2 + 2\z\ 2 + (b + ic)z 2 , \b + ic\ = 2, codimension four (Q4), 

z = —iz + z 2 + (b + ic)z 2 7 generalized Lotka — Volterra (Qf v ), 

z = —iz + z 2 , Hamiltonian triangle, 

where o, b, and c stand for arbitrary real constant. Let 
m ■ = H y (x,y) . = H x (x,y) 

{ ' X M{x,yY V M(x,y)> 

be any of the above systems rewritten in (x, y) coordinates. Here H (x, y) is a 
first integral of system ([T]) with the integrating factor M(x,y). Consider a small 
quadratic perturbations of ([T]): 

/ % • H y (x,y) „ , , . H x (x,y) 

M[x,y) M{x,y) 

where X 2 {x, y, e) and Y 2 {x, y,e) are quadratic polynomials in x and y with coeffi- 
cients depending analytically on the small parameter e. 

Each center of system ([1]) is surrounded by a continuous set of period annuli. 
Compactifying the phase plane K 2 of system ([T]) to the Poincare disc, the boundary 
of the period annulus of the center has two connected components, the center itself 
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and a singular loop which consists orbit (s) and at least one singularity. It is well 
known that the limit cycles of system ([2]) can emerge from 

(a) the center (i.e. the inner boundary), 

(b) the singular loop (i.e., the outer boundary), 

(c) the period annulus. 

Bautin|2] found that at most three limit cycles can appear near a focus or a 
center of any quadratic system. This implies that the cyclicity of the center of 
quadratic system is equal to three under quadratic perturbation. As usual, we use 
the notion of cyclicity for the total number of limit cycles which can emerge from 
a configuration of trajectories (center, period annulus, a singular loop) under a 
perturbation. 

The bifurcation of limit cycles from saddle-loop in perturbations of quadratic 
Hamiltonian systems has been studied in [12] . Moreover, if the loop contains only 
one saddle and under certain genericity conditions, it was proved in [20] that the 
cyclicity of a singular loop can be transferred to the cyclicity of the period annuli. 
However, if the loop contains at least two saddles, this transfer in general is not 
true. For more details, we refer to [18] and references therein. 

The cyclicity of the period annulus of system (JTJ), also known as the (extended) 
infinitesimal 16th Hilbert problem [JJ for n = 2, was investigated by many authors. 
This problem is reduced to counting the number of zeros of the displacement func- 
tion 

(3) d(h, e) = eMt(h) + e 2 M 2 {h) + ■ ■ ■ + M k (h) + ■ ■ ■ , 

where d(h, e) is defined on a section to the flow, which is parameterized by the 
Hamiltonian value h. The number of zeros of the first non-vanishing Melnikov 
function Mk(h) determines the upper bound of the number of limit cycles in (J2J) 
emerging from the periodic annulus of the unperturbed integrable system (fT]) . The 
corresponding Melnikov functions were determined in [14] for quadratic centers. 

The cyclicity of the period annulus for quadratic Hamiltoinian Q$ and Hamil- 
toinian triangle, were completely solved by several authors. See [4j [7l [13l [Til U-3 
|2"31 12~4] and references therein. The generalized Lotka-Volterra Q^ v has been stud- 
ied by Zoladek in [25]. Some results concerned with certain specific case of 
can be found in [3] [6] [TBI HI] etc. However, almost nothing is known about the 
generic reversible case Q§ . Recently, the authors of the paper [9] propose a pro- 
gram for finding the cyclicity of period annuli of quadratic centers of genus one. 
Garu, Manosan and Villadelprat [10] have also some new results in this direction. 

The present paper deals with the cyclicity of period annulus of quadratic codi- 
mension four centers Q4. Using Picard-Fuchs equations and Petrove's method(based 
on the argument principle) [21] , Gavrilov and Iliev[8] proved that the cyclicity of 
period annulus of Q4 is less or equal to eight, see Theorem [2] in Section 2 below. 
Based on their proof, we get the following theorem in this paper. 

Theorem 1. Let system (fT]) be a quadratic codimension four system Q4 rewritten 
in [x, y) coordinates. Then the perturbed quadratic system ^ has at most five limit 
cycles which emerge from the period annulus around the center. Moreover, there 
exists the quadratic polynomials Xi{x,y,z) andYi{x,y,e) such that system ([2]) has 
at least three limit cycles produced by the period annulus of system {TJ . 

The rest of this paper is organized as follows. In Section 2 we briefly sketch the 
proof of Gavrilov- Iliev's Theorem, which is crucial for our analysis. In Section 3 the 
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explicit forms of several related functions are given by revisiting Gavrilov- Iliev's 
proof, and then we get the asymptotic expansions for these functions in Section 4. 
Theorem [T] is proved in Section 5. Finally we give some comments in Section 6. 

2. A SKETCH OF PROOF OF GAVRILOV-lLIEV'S THEOREM 

Gavrilov and Iliev proved the following theorem: 

Theorem 2. [5] The cyclicity of the open period annulus surrounding the center of 
any generic codimension four plane quadratic system is less than or equal to eight. 

We are going to sketch the proof of Theorem [21 

It is well known that the cyclicity of period annulus can be detected in a compact 
region by the number of zeros of the first non-vanishing fc-th order Melnikov func- 
tion Mk(h) in ([3]), which is sometimes called generating function. The generating 
function for Q4 is given in [14, by a complete elliptic integral. After a series of 
changes, the generating function becomes 

(4) 1(h) = Hihl 0fi + /Lt 2 /i,o + M3^o,i + M4(2/-i, + 3KhI- ltl ), 
where 

(5) h,j( h )= [I x*y j dxdy, he (-%—£=), 

JJH(x,y)<h V >>V K / 

with 

(6) H(x,y) = ^(n-l)x 3 -( K -l)x 2 y+^y 3 -y = h, k > 1. 
The integrals U.j in @ satisfy the following Picar-Funchs system 



(7) 



^0,0 


3h 

- 2 ^o,o - 






hjo 


= hl' h0 + 


3^1,1. 




h,i 


_ l r - 

" 3« °'° 


f + 


2(«-l) r/ 

01% 


h,i 


Zh 

- 8 ^0,0 " 








= 3Mi 1)0 


+ 21- 1,1 






K — 1 , 1 , 

= Ii,o + -I- 

K K 


1,0+ 2 l -i,i- 



By using the above system, the authors get 

d_ (I(h)\ hl'(h) - 1(h) G(M 

( ' dh\ h ) h 2 h 2 ' 
where 

(9) G(h) = (^h 2 + /i 3 )^ )0 + + M4[-4W: ii0 + (3nh 2 - 4)I'_ hl }. 
Therefore 



1(h) = h / r 2 <5(£K 

J -2/3 
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and 1(h) has at most as much zeros as G(h). It is proved in [8] that G(h) satisfies 

the following equation 

(10) 

t tu^s d /m a h ((ao + aih 2 + d 2 h* + a 3 h 6 )I^ + (b + b,h 2 + b 2 h*)I{ A ) 
L 2 (h)G(h) = R(h) = , 

where 

d d 2 

(11) L 2 (/i) = 5/c/i - (9nh 2 - 8)— + /i(9k/i 2 - 4) — . 

dh dh z 

Taking the changes 

(12) h = ~JI, Ms) = r 0>0 (h(a)), J 2 (a) = I[ A (h(s)), G(s) = G(h(s)), 
we obtain the equation 

n ^ T r f n sd 2 Id 5\ P 3 (s)Ji(s) + Q 2 ( S )J 2 ( S ) 

where -Ps(s) and Q2(s) are real polynomials of degree at most three and two. Denote 
by dot the differentiation with respect to s. The integrals Ji(s) and J%(s) satisfy 
the following Picard-Fuchs equation 

(14) 6(«-l)(a-/s)' Jl ^ ^ 1_S K-l\fJi 



J 2 ) \\- S S - 1 J \ J 2 

We say that V is a Chebyshev space, provided that each non-zero function in V 
has at most dim(V) — 1 zeros, counted with multiplicity. 

Proposition 3. [5] The following statements hold: 

(i) Suppose the solution space of the homogeneous equation x" + a\(t)x' + 
a 2 (t)x — is a Chebyshev space and let R(t) be an analytic function on 
(a,b) having k zeros (counted with multiplicity). Then every solution x(t) 
of the non-homogeneous equation 

(15) x" + ax(t)x + a 2 (t)x = R(t) 

has at most k + 2 zeros on (a, b). 

(ii) The solution space S associated to the differential operator L 2 (h), defined 
in (jlip . is a Chebyshev space. 

Therefore, we firstly estimate the number of zeros of R(h). Let 

V n = {-P»Ji(s) + Qn-lM*) : Pn,Qn-l € R[s], degP m ,Q m < m} 

Proposition 4. [8] The vector space V n is Chebyshev on the interval (1,k): each 
element has at most dimV"„ — 1 = 2n zeros (counted with multiplicity). 

Proof of Theorem^ It follows from Proposition |4] that P-$(s)Ji + Q 2 (s)J 2 has 6 
zeros in (1, «), and hence R(h) has 6 zeros in (—2/3, — 2/(3(v / k))). Finally one gets 
Theorem [2] from Proposition [3] and (|10JI . □ 

Remark 5. Proposition [7] is proved by using argument principle in the complex 
domain C\(— oo, 1]. The function Ji(s) is a complete elliptic integral of the first 
kind and therefore does not vanish. Let 

. P n { S )Ms) + Qn-l(s)Ms) 

Ji(s) 
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Along the interval (— oo, 1), the increase of the argument of F is bounded by the 
number of zeros of Q n -i(s). Hence we have 

(16) #F(s) < deg P n (s) + #Q„_i(s) + 1. 

in the the complex domain C\(— oo,l], where #F(s) denotes the number of zeros 
ofF(s). 

Remark 6. In the rest of this paper we always suppose k > 1 unless the opposite 
is claimed. For proof's convenience we also suppose that H{x,y), defined in |6j), is 
a first integral of the following system 

dH dH 

(17) X = — = -1 - {K - l)x 2 + ny 2 , y = -— = -2(k - l)x(x - y). 

ay ox 

Hence the annulus Th = {{x , y)\H [x 1 y) — h} has the negative (clockwise) orien- 
tation for h G (— 2/3, — 2/(3-\/k)). The Hamiltonian value h — —2/3 and h — 
— 2/(3y / K) correspond the center (1,1) and the homoclinic loop r_ 2 /(3 v ^) respec- 
tively. 

3. Some paramilitary results 

As in the paper jS], we introduce the variable s € (1, k), defined in (fT2|) . and 
denote by dot the differentiation with respect to s. Taking the changes (fT2"]) and 
(18) 

U«) = I'-i,o( h ( s ))> = (/»(«)), Ms) = Ii, (h(s)), Ms) = Io,Ms)), 
it follows from g]) and Q that 

I(s)=I(h{s)) = (iMi + ^Jx+g^ + fc^Ma)^ 
(19) -6M4y^J3 + (4 + 2s)^ 4 J4 

Suppose that 1(h) is defined as By direct computation we have that, G, 
defined in ©, has the form 
(20) 

Here G(/i) is different from the one defined in ([9]). However if we take 

~ _ 2/x 3 _ _ 2/J.2 2(k - l)(j, 3 
M3 ~ "IT' M2 ~ 3 3^ 

and omit the tildes, then we get G(h), defined in ©. 

For convenience, in what follows we always suppose G(h) is defined in © unless 
the opposite is claimed. 

We note that R(h), defined in (fTTj). has no explicit form in [8 . Following the 
idea in [8] , a direct calculation then yields 

2h(p(h)I>+q(h)I>) 



(21) R(h) 



{9h 2 - 4) 2 (9k/i 2 - 4) 
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where 

(22) p{h) = (9k/i 2 - 4)(o + ai h 2 + a 2 h 4 ), q{h) = b + b x h 2 + b 2 h 4 , 
with 

a = 64^i + 24^ 2 + 8(5K + 3)/i 3 + 32(K-l)/i 4! 

ax = 40(K-9)/xi-162/i 2 -18(ll« + 9)M3 + 24(K-l)(5K-9)/i4, 

a 2 = 18(k + 21)//i + 243k^ 2 + 486k^3 + 54k(k - 1)^4, 

6 = -32((5K-3)/i 2 +3(K-l)// 3 -4(/«-l)/i4), 



(23) 



6i = 72(4(k-1)/ii + (5k 2 + 8k-9)^2 + 3(k-1)(2k + 3)^ 



+4(K-l)(K-3)/i4), 

6 2 = -54(4(«-l)(2K+l)/ii+3«(ll«-9))ti2 + 36«(«-l)/i3 
+12k(k- 1)(2k- 3)/x 4 ). 

Taking the changes (JT3J) , the equation (fl"3)) becomes 
(,4) »W„- «™ 



ds 2 2ds 36/ w 1152(s- k) 2 (s - 1)' 

where 

(25) J-(s) =P 3 (s)Ji(s)+Q2(s)J 2 (s) 
with 

(26) PM = p (-l^,Q 2(s) = q (-lJi 

Here p3(s) and Q 2 (s) are polynomials in s with degp3(s) < 3, degQ 2 (s) < 2. It 
follows from (11411 that 



P'(s) - P 2 (s)Ji(s) + Qi(s)J 2 (s), 

where 

Pi(s) = P 3 (s) — r =a + ais + a 2 s, 

6(s — Kj 

Qi(s) = Q 2 (s) + 6(g _ K){8 _ 1) = A--A), 
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with 

*o = ,),,- "> ~ ; 11 



128(42 + 13k) 16(54+ 13k) 32(27 + 38k + 15k 2 ) 
9^ + " 

128(k- 1)(2k-9) 



3k 



128(-117- 265k + 30k 2 ) 16(-174 + 5k) 16(243 + 121k) 
ai = 27^ Ml + 3^ M2 3^ M3 

64(-1 + k)(-119 + 60k) 

H „ M4, 

9k 



1088(21 + k) 544 1088 1088(-1 + k) 
27k^ k k 9k 



256(-1 + k) 32(-27 + 25k+ 15k 2 ) 16(k - 1)(54 + 31k) 
= 3^ w 3^ ^ 2 3^ ^ 

64(k- 1)(-18 + 5k) 
3^ M4 ' 

64(k- 1)(18 + 77k) 16(— 111 + 137k) 768(k-1) 

A = 27^ Ml 3^ M2 « M3 

64(-1 + k)(-116 + 77k) 

n 

9k 

Solving fii, i = 1,2, 3, 4, from the above equations, we get 

9(162 - 213k - 205k 2 ) 9(23328 + 17604k - 79904k 2 - 59165k 3 ) 

Ml = CKl H OL-2 

P 246400 54454400 

9(27- 18k- 5k 2 ) 9(162 - 4236k + 657k 2 + 2845k 3 ) 

22400(-1 + k) Io 3203200(-1 + k) ^' 

3(54 - 77k + 82k 2 ) 3(3888 + 2502k - 5184k 2 + 11833k 3 ) 

u? = Oil Oil 

P 30800k 3403400k 

3(9 - 7k + 2k 2 ) 3(27 - 709k + 965k 2 - 569k 3 ) 

+ 2800k(k-1) ^° + 200200k(k - 1) l u 

3(-54 + 71k + 205k 2 ) 3(7776 + 5868k - 47598k 2 - 59165k 3 ) 

/iq = ttl Oil 

^ 61600k 13613600k 

3(9- 6k- 5k 2 ) 3(54 - 1412k - 1201k 2 + 2845k 3 ) 

+ 5600k(k-1) ^° + 800800k(k - 1) 1 1 

3(486 - 1017k + 90k 2 + 205k 3 ) 

M4 ~ 246400k(-1 + k) ai 

3(69984 - 162k - 276246k 2 + 44405k 3 + 59165k 4 ) 

H -Cto 

54454400k(k - 1) 

3(k- 3)(-27 + 30k + 5k 2 ) 
p 



22400k(k - 1 
3(486 - 13086k + 16683k 2 + 1050k 3 - 2845k 4 ) 
3203200k(k- l) 2 



8i 
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which yields that 



(27) ao — Po — — noc\ — n 2 a2- 



Remark 7. For proof 's convenience in the rest of this paper we also take a.i,i = 
1,2, and j3i,i = 0, 1, as the new parameters, instead of /J,i, i = 1,2,3,4. 

Without loss of generality suppose fi\ = 1 if f3\ ^ 0. That is to say, j3\ £ {0, 1}. 



Let 



(28) g(s) = ^=P 2 (s) + Q 1 (s)w(s), 

Ji(s) 



where 



(29) u;(«) = P 2 (s) = a 2 s 2 + a lS + a , Q^s) = fts - /3 

Ji(s) 



and a is defiend in ([2T]1. /3 1 £ {0, 1}. 



4. Asymptotic expansions for the related functions 



In this section we are going to give the asymptotic expansions of the related 
functions near the endpoints of their domain of definition. 
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Lemma 8. Ji(s), i = 1,2, ■■■ ,6, have the following asymptotic expansions near 
s = k: 

f 5(s- k) 385(s - k) 2 85085(5- k) 3 
Jl(s) ~ Jl{K) { ~ 36(k - 1) + 5184(k - 1)2 ~ 1679616(k - lj 3 

37182145(s- k) 4 



Ms) = Ji(k) 1 



967458816(k- l) 4 

(s-k) 35(s-k) 2 , 5005(s-k) 3 



36(k - 1) 5184(k - l) 2 1679616(k - l) 3 
1616615(s- k) 4 



967458816(k- l) 4 

T . , A 17(s-k) (1837k -36)(s-k) 2 
!(S) = Jl(K) ( 1 '5i-l| + 5184k(k I) 2 

(5832 - 21276k + 496709k 2 )(s - K f 
1679616k 2 ( k - 1)3 

5(-419904 + 1870128k - 3388824k 2 + 50126789k 3 )(s - k) 4 



+ 



967458816k 3 (k - l) 4 



/ (5k+12)( S -k) (-432 + 1848k + 385k 2 )(s-k) 2 
J4(S) = JlW ( 36k(k-1) + 5184k2(k-1)2 

(69984 - 286416k + 612612k 2 + 85085k 3 )(s - k) 3 

1679616k 3 (k- l) 3 + 

(37182145k 4 + 356948592k 3 - 250327584k 2 + 122332032k - 25194240)(s - k) 4 

967458816k 4 (k - l) 4 

■■■). 

T . , / s-k (36-71k)(s-k) 2 
Ms) = Ji(«) 1 + ^77— TT + 



+ 



+ • 



J 6 (s) = Ji(k) 1 + 



+ 



36(k- 1) 5184k(k- l) 2 

(5832 - 15444k + 14617k 2 )(s - k) 3 
1679616k 2 (k- l) 3 

5(-419904 + 1504656k - 1965816k 2 + 1204387k 3 )(s - k) 4 
967458816k 3 (k - l) 4 

(k - 6)(s - k) (216 - 672k + 71k 2 )(s - k) 2 
36k(k- 1) 5184k 2 (k- l) 2 

(-34992 + 121176k - 185886k 2 + 14617k 3 )(s - k) 3 
1679616k 3 (k - l) 3 

5(2519424 - 10917504k + 18833472k 2 - 19076208k 3 + 1204387k 4 )(s - k) 4 

967458816k 4 (k - l) 4 



+ 
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Proof. Differentiating both sides of system (0, we have 
— T" + I" 

2 o,o ' i o,i — 
Ko + ^i,i 



(30) 



3^0,0 + ^0,1+ 3k 4,i 
g J o,o "I" 2 i»° 4 »i ~ l ~ 4 1,1 



0, 

0, 



3ft/^ 10 + 21^ M 



fC-1 



1 



3// 



T" 

-1,1 



8 J o,o + 4^,1) 

-2/i ll0 . 
1 



Solving I-'j(h) from (|5fj)) . one gets 

^,i 

-1,0 
r" 
-1,1 

V Io,i J 



(31) 



(9ft 2 -4)(9k/i 2 -4) 



M 



^-i,i- 



f \ 

Ii.i 

v ^1,1 y 



M = 



where 

/ -3ft(9«;ft 2 -4) 12(k- l)ft 

-3h(9nh 2 - 4) 3/j(9k/i 2 - 4) 

8(k-1) -8(/c-1) 

-12(K-l)ft 12(«-l)ft 

2(9«ft 2 - 4) -2(9k/i 2 - 4) 

\ 2(9 K ft 2 -4) -18(K-l)ft 2 

It follows ([21]), (H2) and (0 that 

/ Ji \ 






6/vft(9/i 2 - 4) 2 K (9ft 2 - 4) 

4(9ft 2 - 4) -3«ft(9ft 2 - 4) 



(32) 



with 



6(-l)(,-«) 5 



J 2 



M* 



■h 

V ^ / 



M* = 



1 - s 




K - 1 


1 - s 




s- 1 


\/k{k - 


1) 




V~s 






1 — K 




K- 1 


\/k(s - 


1) 


v^(s - 1) 








\/k(s - 


1) 


(k - 1)s 






-2(s-k) - 
2(s - k) 






^k(s - k) 



K — S 
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Since h = —2/3 corresponds the center (1,1) of Hamiltonian system (fl7|). we have 
I, j (-2/3) = 0, i = 0,±1, j = 0,±1. Substituting X iJ (-2/3) = into Q yields 
/l M (-2/3) - /l li0 (-2/3) - J( i0 (-2/3) - = ^(-2/3) - ^, (-2/3). 

Therefore Ji(k) = J 2 (k) = • • • = Jq(k). 

Since s = K corresponds to the center of Hamiltonian system (|17l) . Ji(s) is 
analytic at s = k. Taking Ji(s) = J2fLo Ci ,j( s ~ with c^.o = Ji(rt) into (l3"2"j) . we 
get the expansions. □ 

Corollary 9. TTie following assertions hold: 

(i) J(s) = if and only if A*i = A*2 = fJ.3 = A*4 = 0. 

(ii) There exists i — 1,2,3,4, such that I(s) has at most three zeros in 
(1,k). 

Proof. It follows from (fl"9)) and Lemma [8] that I(s) has the following asymptotic 
expansion at s = k: 

(33) I(s) =v 1 (s-k) + v 2 {s - k) 2 + v 3 (s - k) 3 + v 4 (s — k) 4 H , 

where z/j = d l I (k) / ds z , i = 1,2,3,4, with 

2a* 1 A*2 A*3 . 2(«- 1)a*4 



v 2 



9k 3k 3k 3k 

(13k- 18)a*i (17k- 18)a* 2 (17k - 12)a<3 (13k + 18)aJ4 
162k 2 (k- 1) + 108k 2 (k- 1) + 108k 2 (k- 1) + 54k 2 ' 
_ (1944 - 4068k + 1739k 2 )a*i (1944 - 3780k + 1801k 2 )ai 2 
11664k 3 (k- l) 2 (7776k 3 (k- l) 2 

(1296 - 2712k + 180lK 2 )Ai 3 (-1944 + 1260k + 1739k 2 )a(4 



(34) 

v 4 



7776k 3 (k- 1) 2 3888k 3 (k-1) 2 
5(-104976 + 324648k - 338526k 2 + 101837k 3 )a*i 
1259712k 4 (k- l) 3 
5(-104976 + 309096k - 301374k 2 + 96253^2 
+ 839808k 4 (k- l) 3 

5(-69984 + 216432k - 225684k 2 + 96253K 3 )Ai 3 
+ 839808k 4 (k - l) 3 

5(104976 - 180792k + 7398k 2 + 101837k 3 )aj 4 
+ 419904k 4 (k - l) 3 ' 

System (|34[) is a linear system of equations in the variables A*i, i = 1,2,3,4. The 
determinant of matrix of coefficients of (|34l) is equal to 125/(472392k 8 (k — l) 4 ) > 
for k > 1. As shown by Cramer's rule, system Q34I) has a unique solution. Therefore 
V\ = v% = v% = v± = if and only if a*i = A*2 = A*3 = A*4 = 0, which yields that 
I(s) = if and only if /ij = 0, i — 1, 2, 3, 4. This proves the assertion (i). 

Since system (|34[) has a unique solution, we can choose i>i as the independent 
parameters, instead of A*i, i = 1, 2, 3, 4. Denote by 7(s, ^1, 1^2, ^3, ^4) = -f (s). With- 
out loss of generality suppose 1^4 > 0. To get more zeros of Z(s), we choose t'i 
and Si G (1,k),i = 4,3,2,1, such that I(si, 0, 0, 0, 1/4) > 0, 1(53,0,0,^3,^4) < 
0, 7(s 2 , 0,^2,^3,^4) > 0, 7(si, 1/1,1*2, z/3, ^4) < and < |^i| < \v 2 \ < |^| < M, 
1 < S4 < S3 < S2 < s% < k. It is easy to show that I(s), defined as the above, has 
at least three zeros which tend to k. □ 
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Lemma 10. Ji(s), i = 1,2, ••• ,6, have the following asymptotic expansions near 
s = 1; 

ln(.s-l) 5(s-l)ln(s-l) 
Jl(5) = -V^ = T +12_ 72(«- 1)3/2 

3 (s-l)ln(s-l) 



J 2 ( S ) = 



12(k- l) 3 / 2 



where C\% is a real constant. 

Proof. Since the period annuli of the vector field (fl~7j) begin at the center (1,1) 
and terminate at a homoclinic loop r_ 2 /(3 v /^) = {{x, y)\H (x, y) = —2/(3-^/^)}, it 
follows from [35] that Ii j(h), i > 0, have the asymptotic expansions of the form 



ft 



hi 



as ft — > —2/(3s/k), which implies that Ji(s), i = 1,2, have the asymptotic expan- 
sions of the form 

(35) Jj(s) = Ci,i ln(s - 1) + Cf, 2 + Ci >3 (s - 1) ln(s - 1) + • • ■ 

as s — > 1. On the other hand it is well known (see for instance JT71 [23] or the 
appendix of [24]) that 



x l y J 



20^1 



1 



= <^ 2 a /kJ(«-1) 



0. 



if i = 0, 
if i = 1. 



(x, ! /)=(0,l/ x /S) 

A simple calculation shows that c^i = — dij t i. Taking (|35[) with c^i = — <2ij,i into 
(|14[) . we obtain the asymptotic expansions near s = 1 for Ji(s), i = 1, 2. □ 



By Lemma [8] and Lemma [101 a straight calculation shows that the following two 
lemmas hold. 



Lemma 11. The following expressions hold: 

1 ... s 25 



w(k) = 1, «/(/s) = 



6(k-1) 



it) (ft) = — 



216(k- l) 2 ' 



u/"(k) = 



775 



3888(k- 1) 



fl («) = 0, T{k) = F{k) = 0, T"(k) = Ji(K)g'(K). 

Since G(s) is analytic at s = k, it follows from ([24|) that .F(k) = F'{k) = 0. This 
is verified by Lemma 1111 

Lemma 12. The the following expansions holds as s — > 1: 
w(s) = 



6 



lis) = 



ln(s - 1) 
' ft(l) - 



6Qi(l) 
ln(s - 1) 



6(s- 1) 
ln(s - 1) 



if Qi(l) ^0, 
if /3 = A = 0, 

if A, = A = 1, 
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where Pi £ {0, 1}. 

5. Proof of Theorem [T] 
First of all we note that the following proposition holds: 

Proposition 13. Denote by #T{s) the number of zeros ofJ-(s). Then we have 

(36) #F{s) < #P{s) = #g{s), #!(«) < #G{s) < #F{s) + 2. 

Proof. It follows from Lemma [11] that J-(n) = J-'(k) = This yields the first 
inequality of (f3l)|). The second inequality is obtained by Proposition [3] □ 

Noting that the cyclicity of period annulus is determined by #/(s), we will prove 
Theorem [1] by estimating the number of zeros of g(s). 

Since g(n) — 0, we get that g(s) has at most three zeros in (1,k) by argument 
principle, see Remark [5] This implies that #/(s) < #G(s) < 5. However, to get 
more information about the number of zeros of g(s) (hence I(s)), we prefer to prove 
Theorem [T] by the following theorem, see the comments in the next section, and the 
note after the statement of this theorem. 

Theorem 14. Let s £ (1, «) and pi £ {0, 1}. 

(a) If pi = 1, /3q £ (—oo, (23k — 54)/31] U [1, +oo), then g(s) has at most two 
zeros in (1, k). 

(b) If f}\ = 1, flo £ ((23k — 54)/31, 1), then g(s) has at most three zeros. 

(c) If Pi = 0, Po 7^ 0; then g(s) has at most two zeros in the same interval. 

(d) If Pi =0, Pq — 0, then g(s) has at most one zeros. 

If pi = 1, Pq £ [l,+oo) (resp. Po £ (— oo,l)), then it can be proved that 
g(s) has at most two (resp. three) zeros in (1,k) by argument principle, see the 
proof of Proposition [4^8] , or Remark [5j However it seems that we can not prove 
by argument principle that g(s) has at most two zeros in (1,k) if Pi = 1, Pq £ 
(-oo, (23k- 54)/31]. 

Firstly we study the geometric properties of w(s) = ^(s)/ Ji(s). 

Lemma 15. The function w(s) is monotonically increasing and concave in the 
interval (1, k), i.e., w'(s) > 0, w"(s) < and < w(s) < 1. 

Proof. It follows from (jT4j) that w(s) satisfies 

(37) 6(s - l)(s - k)w' = 1 - s + 2(s - l)w - (« - l)w 2 = U{s, w). 

Note that U(s,w), the right hand of (|37|) . is a quadratic polynomial of w. Since 
4(s- 1) 2 + 4(k- 1)(1 -s) = i(s-l)(s-K) < for s £ (1,k) and -(k- 1) < 0, we 
have U(s, w) < 0. This yields w'(s) > for s G (1, «). The inequality < w(s) < 1 
follows from Lemma [Tl] and Lemma [12] 

Differentiating both sides of ([37]) . one gets 



6( S - 1)( S - «)» = (-6(2. - k - 1) + -g— j 6(g _ 1)(g _ K) + — g^— , 
which implies that 

d 2 w _ Vi(s,w)V 2 (s,w) 



(38) 



ds 2 18(s- l) 2 (s- k) ; 
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with 

Vi{s,w) = (k - 1)w- (s- 1), V 2 {s,w) = («- l)w 2 + (4s - 3k - l)w - 2(s - 1). 

By Remark[5](or the proof of Proposition S[8 ), we conclude that (1 — s) Ji(s) + (k — 
1)>^2(s) has at most two zeros in C\(— 00, 1]. Since (1 — s)J±(s) + (k — l)J 2 (s) has 
a zero at s = k, (1 — s) Ji(s) + (k — 1) J2(s) has at most one zero in (1, ft). Noting 
Ji(s) ^ 0, we know that the function 

m ( s ) = V 1 (s,w(s)) = (l-s) + (n-l)^f\ 

Ji(s) 

has at most one zero in (1, k). By direct computation 771 (k) = r\\ (1) = 0, lims^i 77^ (s) 
+00, Tj[(k) = —5/6, which implies that the number of zeros of rji(s) is even. There- 
fore 771 (s) has no zero in (1, k), which shows that V\(s, w(s)) > in (1, k). 
Let 

(39) V = {(s,w)\l < s < k, < w < 1}. 

Now we study the two independent variables function Va(s, if), defined in T>. The 
equations dV 2 /ds = dV 2 /dw = has a unique solution at (s, w) = ((/t+l)/2, 1/2) £ 
T>. Therefore, V 2 (s,w) has a maximum point and a minimum point at either 
(s,w) = ((k + l)/2, 1/2) e T>, or the point in the boundary of T>. Since 

V 2 {s, 0) = -2(s - 1) < 0, V 2 {s, 1) = 2(s - k) < 0, ^a(l, i«) = (/s- l)w(w - 3) < 0, 

V 2 (k,w) = ( K -1)(w-1)(w + 2)<0, V 2 (Z±±,~) =-^^<0, 

we have V 2 (s, w) < for (s, iu) e V. 

Let 7y 2 (s) = V2(s,w(s)). Then 772(1) = tj 2 (k) = 0, lim s _,.i T)' 2 (s) = -oo,?7 2 (k) = 
5/2 > 0, which implies that V 2 (s, w{s)) < in (1, n). It follows from (J38J) iu"(s) < 
for s € (1, k). □ 

Assume that i()(xi,x 2 , - • • ,x n ) and (f>(xi, x 2 , ■ ■ ■ ,x n ) are two polynomials in 
x\,x 2l --- ,x n . Eliminating of the variable xt from the equations ip(xi, x 2l ■■ ■ ,x n ) 
4>(xi,x 2 , • • • , x n ) = 0, we get the resultant of ip(x±, x 2 , ■ ■ ■ , x n ) and 4>{x\,x 2 , • • • , x n) . 
denoted by Resultant (ip, <j>,Xi). 

Lemma 16. w"'(s) > 0, s £ (1,k). 
Proof. It follows from (38|) and {37]) that 

$(s, 10) 



u/"(s) 



108(s-1) 3 (s-k) 3 ' 
with 

$(s,w) = -(s-1) 2 (20s + k-21) + 2(s-1)(15-k + 6k 2 -29s-11ks + 20s 2 )w 
-2(k - 1)(1 + 18k - 19s) (k - s)w 2 + 6(k - 1) 2 (1 + 3k - 4s)w 3 
-3(k-1) 3 ti; 4 . 

Let I? be the closed rectangle, defined in (|39[) . The maximum and minimum for 
$(s, u>) in V necessary occurs either on the boundary of V, or the points inside V 



PERTURBATION OF Q 4 



15 



whose coordinates satisfies equations 

40) <J> S = = 0, 3> w = — = 0. 

as aw 

4> s and § w , denned in (|40p , are two polynomials of s with the polynomial coefficients 
depending on w and k. Their resultant is 

Resultant ($ a ,$ w> s) = -8000(k- 1) 6 (u; - l) 2 w 2 (2w - l)x(w) 

with 

X(w) = -160425 + 316012™ - 314956u; 2 - 2112w 3 + 1056w 4 . 

It is nice for our study that Resultants ($ s , <$> w , s) does not depend on k. By Sturm 
Theorem x( w ) h as n0 rea l zero in (0, 1). Taking w = 1/2 into the first equation of 
system (|40|) . we know that (s,w) = ((fc + l)/2, 1/2) is a solution of system (|40|) in 
T>. Direct computation yields that if (s,w) G T>, then 

\ 2 '2/ 16 V ; 

$(l,tu) = -3(k- 1) 3 w 2 (12-6w + w 2 ) < 0, 

$(k,w) = -3(k-1) 3 0-1) 2 (7 + 4w + w 2 ) < 0, 

$(s,0) = -(s-1) 2 (20s + k-21) <0, 

$(s, 1) = -(s - k) 2 (-1 + 21k - 20s) < 0, 

which imply that the maximum and minimum for $(s,w) are non-positive. There- 
fore $(s, w) < for (s, w) e P. 

Assume that there exists the internal point (s*,w*) of V such that $(s*,u>*) = 
0. Since $(s, iu) < 0, (s*,u>*) must be a maximum point of $(s, iu) inside P. 
However we have shown that the maximum and minimum for <£>(s,ui) inside D 
necessary occurs at (s, w) = ((k + l)/2, 1/2) and $((« + l)/2, 1/2) < 0. This yields 
contradiction. Hence <S>(s,w) < for (s,w) € {(s,u;)|l < s < /c, < w < 1}, which 
implies that w"'(s) > for s e (1, k). □ 

Proposition 17. Let /3i = 1 and s € (1, as). TTie following statements hold. 

(i) #A) G (-00,(54 -23k)/31], tfien g"'{s) > 0. 

(ii) 7/ Pq € ((54 — 23k)/31, 1), f/ien g"'(s) has exactly one zero. 

(iii) J/0„ € [l,+oo), then g"'(s) < 0. 

Proof. By direct computation we have 

g"'{s) = 3w"(s) + (s- poW"{s) = 3w"'(s)6(s), 

where 

Q(s) _ «"(-) , s-/3 



to'"(s) 3 ' 
Therefore, 

d0(«) _ 4(u/"(s)) 2 -3ii/'(s)u/ 4 >(s) _ e 1 (s,w,K)e 2 (s,w,K) 



(41) 



ds 3(w"'(s)) 2 17496(s - l) 6 (s - k) 6 w"' 2 (s) 
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with 

9 1 (s,«;,k) = 2(s-1) 2 (-9 + 4k + 5s)-(s- 1)(36- 67k + 51k 2 -5s-35ks 
+20s 2 )w + (k - 1)(18 - 41k + 18k 2 + 5s + 5ks - 5s 2 )w 2 
± 6 a ( S )+9 1 ( S )w + d 2 (s)w 2 , 

@ 2 (s,w,k) = (s-1) 2 (9 + 7k-16s)-(s-1)(-9-62k + 39k 2 + 80s + 16ks 
-32s 2 )w - (k - l)(-9 + 55k + 18k 2 - 37s - 91ks + 64s 2 )w 2 

+9(k- 1) 2 (1 + k-2s)w 3 . 

A calculation shows that #2(1) = ^(k) = 18(k— l) 3 > 0. Since 6*2(3) is a quadratic 
polynomial in s with \im s ^± 00 02(s) = —00, we have #2(5) > for s £ [1,k]. On 
the other hand, 

(0i(s)) 2 - 40 o (s)0 2 (s) = 25(s - k) 2 (s - 1) 2 (81 - 146k + 81k 2 - 16s - 16ks + 16s 2 ). 

It follows from (16 + 16k) 2 - 4(81 - 146k + 81k 2 ) • 16 = -4928(k - l) 2 < that 
81-146k + 81k 2 -16s-16ks + 16s 2 > 0, which implies (6 1 (s)) 2 - 40 o (s)0 2 (s) > 0. 
Rewrite 0i(s,iu, k) as the form 

01 (a, 10, k) = 9 2 (s)(w - W + (s))(w - W~(a)), 

where 

± _ -9,(8) ± v/(0i(s)) 2 - 40o (s)0 2 (s) 



202 (a) 
This gives that 



3(k- 1) ' w 6(k- 1) 

1+, and 

(s-k) 115(s-k) 2 



W + {s) = 1 
W~{s) = 1 



6(«-l) 324(k-1) 2 
8(s - k) 



3(k-1) 

as s — > k~. Therefore it follows from Lemma [TT] and Lemma [T2l that W~(s) < 
W + (s) < w(s) — J 2 (s)/Ji(s) as either s — > 1 + , or s — > n~. 

Since W + (s) < w(s) as either s — > 1 + , or s — > k~ , w = w(s) and w — W + (s) 
have at least two intersection points if there exists. Suppose that w = w(s) inter- 
sects w = W + (s) at si, S2, then w(si) — W + (si) = 0. By the Mean Value Theorem 
there exist s* such that w'(s*) - (W+)'(s*) = 0. Noting (VK ± )'(s) = -G u /9iw,ti 
follows from (|37|) that s — s* satisfies the following equations 



(42) 6i(s,w,k)=0, 9i(s,ii),K) =@ lw U(s,w) +6(s -l)(s- k)Q 1s , 

where Oi s = d<di(s,w,k)/ds. Noting that we have shown 02(s) > for s £ (1,k), 
one gets 

Resultant (61,61,-0;) = -35083125(k- l) 5 (/c - s) 5 (s - l) 5 2 (s) ^ 0, s e (1,k), 
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which implies that two equations in (|42|) have no common zero. Therefore there is 
no s* G (1, re) such that w(s*) = W+(s*). This yields W-(s) < W+{s) < w(s) for 
s £ (1,k). Finally we obtain that Q 1 (s,w(s),k) = 28 2 (s)(w(s) - W + (s))(w(s) - 
W~(s)) > for s £ (1,/t). 

Now we consider 2 (s, w, k). Let 

V' = {(s, k)|1 < s < K, 1 < k < c}, c is a real large constant, |c| S> 1}, 

be a triangle in the res-plane and fix w as a real constant with < w < 1. The max- 
imum and minimum for 2 (s, it), re) in I?' necessary occurs either on the boundary 
of V , or the points inside V whose coordinates satisfies equations 

d<d 2 (s,w,K) _ dQ 2 (s,w,K) _ 
ds dk 
Direct computation shows that 

Resultant(0 2s , 2k , w) = -6705(re - l) 2 (s - re) 2 (s - l) 2 -f(s, re) 

with 

7(3, k) = 362313 + 701586k- 1012697k 2 + 421884k 3 - 1012697k 4 + 701586re 5 
+362313k 6 - 8(1 + k)(359433 - 174116k - 174026k 2 - 174116fc 3 
+359433k 4 )s + 8(991241 + 22492k - 1044426k 2 + 22492/fc 3 
+991241k 4 )s 2 - 16384(k + 1)(649 - 978k + 649k 2 )s 3 

+8192(809 - 658k + 809k 2 )s 4 - 1572864(1 + k)s 5 + 524288s 6 . 

Since Resultant (7,5, 7 K , s) = c*(k— l) 25 with c* < 0, the maximum and the minimum 
for 7(s, k) in V occurs on the boundary of V . 7(s, k) is a polynomial of k with 
degree 6 and the coefficient of k 6 is 362313, which implies 7(s, c) > as c is sufficient 
large enough. Noting 7(1, k) = 7(k, k) = 362313(k — l) 6 , 7(3, k) has its minimum 
value 7(1, 1) = at (s, k) = (1, 1). This yields 7(3, k) > for (s, re) G P'\(l, 1). 

Therefore Resultant(6 2;s , 9 2re , to) 7^ for (s,k) G 2?'\({s = l}U{s = re}), which 
implies that the maximum and minimum for 2 (s, w, k) in T>' necessary occurs on 
the boundary of T>' . If < w < 1 and (s, k) 7^ (1, 1), then 

62(1,10, «) = 9(k-1) 3 u; 2 (w-2) < 0, Q 2 (k,w,k) = -9(re- lf{w - l) 2 (w + 1) < 0. 

Noting that 2 (1, w, re) is a polynomial in k and the coefficient of the highest order 
term k 3 is 9w 2 (w~ 2) < 0, we have 2 (s, u>, c) < 0, provided that c is sufficient large 
enough and < w < 1. Summing the above discussions and noting 2 (1, w, 1) = 0, 
one gets 2 (s, to, k) has its maximum value zero at (s, k) = (1, 1) in D' . Since we 
always suppose that k > 1 in this paper, 2 (s, w(s), re) < for s £ (1, re). 

It follows from fill) that 0'(s) < 0. This yields that g"'(s) has at most one zero 
in (1, re). On the other hand, Lemma [TT1 and Lemma [12] give 

r ,„, , / +°°' if ^ < X ' 25(23re + 31/j -54) 

limg s) = < .g fc = — — -3 , 

{ -00, if fa > 1, 3888(re- l) 3 

which implies the assertions of this proposition. □ 



Corollary 18. Let s £ (1,k) and fix = 1. 
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(a) If Po £ (— oo, (23k — 54)/31] U [1, +00), then g(s) has at most one inflection 
points. 

(b) If /3q € ((23k — 54)/31, 1), then g(s) has at most two inflection points. 

Proof. Note that the zero of g"'(s) is the maximum or minimum point of g"(s) and 
g(n) = 0. The assertions of this proposition follows from Proposition [T71 □ 

Proof of Theorem \14\ Firstly we note that g(n) — 0. 

If /?i = 1, then the statement (a) and (b) follows from Corollary [T51 

If (3i = 0,ft ^ 0, then Lemma Ql shows that g'"(s) = /3 w"'(s) ^ 0, which 

implies that g(s) has at most one inflection point. This yields the statement (c). 
If Po — Pi = 0, then g(s) — P2(s) = (s — n){ai + na.2 + a^s). The assertion (d) 

follows. □ 

In the end of this section, we prove Theorem [TJ 

Proof of Theorem^ Proposition [TBI and Theorem [T4l show that I(s) has at most 
five zeros in (1,/c). This implies that the perturbed system |(3J) has at most five 
limit cycles which emerge from the period annulus around the center. The second 
assertion of Theorem [JJ follows from Corollary [5] □ 

6. Comments 

Zoladek conjectured that the exact upper bound of the cyclicity of the period 
annulus for Q4 is three [14j [25] . Unfortunately we can not prove Zoladek's conjecture 
in this paper. 

As we mentioned before, the argument principle gives a shorter proof of Theorem 
[1] However it seems clear to us that it does not allow to go further in Zoladek's 
conjecture. Our approach is perhaps more involved from the computational point 
of view, but we think that it may provide a way to attack the problem in the future 
paper. For instance, we can get the following results from Lemma II 5[ Lemma [TBI 
and Corollary IT81 

1. If pi = 1, fa < 1, Wo) < ( res P- A) > k, P 2 (/? ) > 0), then {P 2 {s)/Q 1 {s))"+ 
w"(s) = P2(Po)/(s — fio) 3 + w"(s) < 0. This implies that g(s) has at most one zero 
in (1,/c). Therefore I(s) has three zeros in the same interval. 

2. If /3i = 1, fa < 1, < P 2 (Po) < 25(1-/3 )/(432(k-1) 2 ), then (P 2 (s)/Q 1 (s))"+ 
w"(s) = P 2 (Po)/{s - Po) 3 + w"{s) < 0. This yields J(s) has three zeros in (1, «). 

3. Ifp l = l, po G (-oo,(23k-54)/31]U[1,+oo), P 2 (l)g'(K) > 0, then g(s) has 
at most one zero. Hence I(s) has at most three zeros. 

Here we just list the partial results we have proved. We wish that the above 
results will be helpful for proving Zoladek's conjecture in the future paper. 
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